Generation of optical combs in a whispering gallery mode resonator from a 

bichromatic pump 
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An optical comb is shown to arise from a whispering gallery mode resonator pumped by two optical 
frequencies. Two externally excited modes couple due to Kerr nonlinearity to initially empty modes 
and give rise to new frequency components. This thresholdles process is much more efhcient than the 
previously reported single-pump four-wave mixing. As a result, a few milliwatt pump is sufBcient 
to generate strong secondary fields, that efficiently generate higher-order frequency components and 
so on, in a cascade process leading to an optical comb. 

PACS numbers: 42.60.Da, 42.65.Ky, 42.65.Hw 



Optical frequency combs 1_, _2] found their applications 
as optical clocks and frequency standards of very high ac- 
curacy 3i| , including those for astronomy [4| and molecu- 
lar spectroscopy [5, 6] . They have also been proposed for 
use in quantum information processing [7|, |8| . Four- wave 
mixing in a resonator with Kerr nonlinearity |9J, llO| has 
been recently demonstrated [ll|,[l3] as one of the methods 
to generate optical combs. In this case the comb arises 
from direct or cascaded frequency conversion processes, 
which requires very high effective nonlinearities. The lat- 
ter can be achieved by using the whispering gallery mode 
(WGM) resonators, known for their high Q factors. 

In monochromatically pumped Kerr-based optical 
comb sources two pump photons are converted into a 
quantum-correlated photon pair. This process, whose en- 
ergy diagram is shown in Fig.[T]^a), populates the pair of 
modes adjacent to the pump as well as those two, three, 
etc., of the resonator free spectral ranges (FSR) away 
from the pump. To generate such a frequency comb in a 
WGM resonator, the pump has to exceed a certain power 
threshold Pth, inversely proportional to Q^ [l^. The fac- 
tor Q^ also determines the phase drift of the comb com- 
ponets, via the Schawlow-Townes formula [l3|. On the 
other hand, the frequency span of the comb is limited by 
dispersion and therefore scales as the resonator linewidth, 
i.e. as 1/Q. This leads to a trade-off between the comb 
pump-efficiency and stability, and its span. 

In this Letter, we report a different configuration of 
a four-wave mixing comb using two pumps coupled to 
two WGMs one or several FSRs apart. Instead of cou- 
pling two vacuum and two (degenerate) pump fields, the 
leading-order interaction couples three pump fields and 
one vacuum field. For example, on the energy diagram 
in Fig.[ljb) two pump photons with frequency wi are ab- 
sorbed, and one pump photon with frequency lj2 > wi 
is emitted, along with the photon with u;_ = 2wi — uj2. 
A symmetric process will lead to generation of a pho- 
ton with Lj^ = 2lj2 — LOi- Since the pump field is much 
stronger than the vacuum field, the non-degenerate four- 
wave mixing is expected to be much more efficient than 
the degenerate one. Furthermore, as we will see below, 
the non-degenerate process is thresholdless. It occurs at 



any pump power, which allows one to avoid the undesir- 
able high-power effects such as thermorefractive oscilla- 
tions [13l| and stimulated Raman scattering. 

High efficiency and absence of a threshold lead to a cas- 
cade process in which the new frequency components gen- 
erate further components, and a frequency comb arises. 
The spacing of this comb is equal to W2 — ^i , which can 
be a multiple of the FSR. The relative frequency drift of 
this comb components is determined by the drift of the 
lasers beat note, and is independent of Q, which resolves 
the stability - span dilemma. Phase-locking the lasers 
beat note to a high-stability reference oscillator one can 
achieve high comb stability even with low-Q resonators 
that are suitable for broader-range combs. 
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FIG. 1: Energy diagrams of previously (a) and presently re- 
alized (b) four-wave mixing processes generating a frequency 
comb. The experimental setup (c). 

Our experimental setup diagram is shown in Fig.[IJc). 
We used a fluorite resonator with 13.56 GHz FSR. Fluo- 
rite has nearly linear dispersion at the wavelength of 1.5 
/im, so the resonator spectrum around this wavelength is 
highly equidistant. Light was coupled in and out of the 
resonator using two optical fibers polished at the optimal 
coupling angle [ij]. The light from the input fiber re- 
flected from the resonator's rim was collected by a photo 



detector to observe the spectrum of the resonator. 

The input fiber combined hght from two lasers centered 
at 1560 nm. Both laser's frequencies were simultaneously 
scanned around the selected WGMs of the same family. 
This was achieved by fine-tuning each frequency until 
the selected resonances observed by the photo detector 
overlap. The resonator quality factor Q — 10® (loaded) 
was made relatively low to increase the linewidth and 
to make the overlap easier to achieve. An optical spec- 
trum analyzer (OSA) is connected to the output fiber for 
continuous data acquisition. The OSA was set to retain 
the peak power values, therefore a trace recorded for a 
sufficiently long period of time represented the situation 
when both lasers are fully coupled to the WGMs. 
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FIG. 2: A typical spectrum of a bichromatically-pumped flu- 
orite WGM resonator. The two highest peaks correspond to 
the two pumps, each of 2.9 mW at the input. 

A typical comb generated by two pumps of -Pi = -P2 = 
2.9 mW set one FSR apart is shown in Fig. [21 We did 
not see any pump power threshold required for the on- 
set of the nonlinear oscillations. This measurement was 
repeated for the laser frequencies set two, three and ten 
FSRs apart. We had to stop at ten FSR because the 
lasers could not be tuned further. No significant differ- 
ence between these measurements and the first one was 
observed. A comb spectrum corresponding to the ten 
FSR separation and the pump power Pi = P2 = 4.75 
mW is shown in Fig. [31 The output pump power in 
Figs. [21 and [3l is much less than the input powers be- 
cause the output coupling was made weak so as not to 
overload the resonator. 

Our experiment can be described within the general 
theoretical model of Kerr-couphng of the WGMs p^ . 
The difference in the present treatment is that the par- 
tial degeneracy of the Hamiltonian in [ij] is removed, 
and the leading-order process is a four-mode interaction 
rather than three-mode. The Hamiltonian H = Hq + V 
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FIG. 3: The spectrum for Pi : 
(135.6 GHz) pump separation. 



P2 = 4.75 mW and 10 FSR 



describing our system includes the following terms: 

V ^ -h^ ■.(a + b + c + d + h.c.)'^:, (1) 



-h^: {a + b + c 



where c, a, b, and d are photon annihilation operators 
at frequencies ujc < oja < uii, < uJd, respectively. The 
coupling constant g is found in [13l |: 
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(2) 



where 712 is the Kerr nonlinearity of the resonator mate- 
rial, no is its refraction index and V is the mode volume. 
Notice that the fields represented by the operators a 
and b are strong pump fields, while the fields represented 
by c and d are initially vacuum fields. In Hamiltonian 
([T|) we will only retain the leading terms with at least 
three operators representing strong fields. Following the 
steps of '13' and introducing the Langevin formalism, we 
arrive at the following stationary state equations: 



(3) 



where the loss rate 7 and optical frequency ui are assumed 
to be the same for all modes. The Langevin forces 



7a = 


= 6ig{a''aa + 2ab% + 2a''bc + bbd^) + fa, 
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(4) 



represent the external pump coupled into the WGMs. We 
have assumed critical input coupling, so that the coupling 
rate 7 in (|3]) and the loss rate 7 in ([3]) are equal. As we 
have mentioned, the output coupling rate 70^4 is much 



less than 7, so we have neglected the terms with ^out in 
Eqs. dS]). The output coupling rate determines the rela- 
tion between the average of the intracavity field operator, 
e.g. (a), and the output power supplied to the OS A: 



Pr* = {afloutnu/2. 



(5) 



In the following we will assume equal pump powers: 
Pa = Pb = P- The solutions to the first pair of Eqs. ^ 
retaining the leading-order self- and cross-phase modula- 
tion terms are 
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(6) 
(7) 



where F{P/Pq) is the real root of the cubic equation 
9^3 _|_ ^ _ P/Pq = 0. The scaling power Pq is expressed 
via the single-pump hyperparametric oscillations thresh- 
old power Pth, introduced in [13]. Assuming A ~ 1.56 /im, 
no = 1.44, n2 = 3.2 • IQ~^^ cmVW, Q « 10^ V « lO^* 
cm'^, we find for our resonator Pq ~ 50 mW. 

Remarkably, the low-power limit of the function 
F{P/Po) is P/Pq, which is also the solution of the first 
pair of Eqs. ([3]) with all nonlinear terms neglected. At 
higher powers F(P/ Pq) grows slower than P/Pq, which is 
indicative of the "pumping inefficiency" [T^ due to self- 
and cross- phase modulation of the pump fields. 

From the last pair of Eqs. ([3]) we find, using Eqs. ([6]) 
and (O, the intracavity field operators mean values and 
the output powers for the new optical fields. Notice that 
the following solutions have no pump power threshold: 

(c) = id) = ^-l{af = j2.F'/'{P/Po), (8) 
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(9) 



The output pump and signal powers delivered to the 
OS A are given by Eqs. ([5]) and ([5]), respectively. This 
allows us to determine F{P/Pq) at e.g. P — 0.75 mW, 
and therefore to find the experimental value of Pq. We 
have carried out this measurement for one, two, three and 
ten FSR pump separation and foimd Pq = 10 ± 1 mW. 
This value is a factor of five below our theoretical esti- 
mate. The discrepancy is likely due to large uncertainty 
in the shape of the hand-polished resonator, and the con- 
sequent order-of-magnitude uncertainty in the mode vol- 
ume V. Also, the intrinsic Q of the resonator could be 
larger that the measured loaded value. 

To compare the theoretical prediction ^ with the ex- 
periment we have measured the ratio P^^\P) / P'-^^ {P ~ 
0.75 mW) and plotted it as a function of the input pump 
power P, see Fig.|4l Also shown are the theoretical curve 
(P|) and the cubic power curve, which would be the so- 
lution in the absence of the "pumping inefficiency" . We 
see that the role of this inefficiency is significant. 
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FIG. 4: Normalized power of the first-order frequency com- 
ponents vs. the pump power: the data and the theory curves. 

We have already mentioned that the bichromatically 
pumped process is thresholdless and therefore can be ef- 
ficiently cascaded, with each excited mode playing the 
role of a new pump. Carrying out the above analysis re- 
cursively for the higher-order comb components, we find 
the following relation between their powers: 



p(n)/p(l)=^,^(^p(2)/p(l) 



(10) 



The result PH)) describes the comb shape in the nor- 
malized form that can be easily verified in experiment. 
The factor ^„ has been introduced to account for mul- 
tiple leading-order channels available for generation of 
sufficiently high-order frequency components, see Fig. [5l 
Each channel shown in Fig. [S] annihilates a photon pair 
in some blue-shifted mode m and couples it to a yet un- 
populated blue-shifted mode n > m and to already pop- 
ulated red-shifted mode n — 2m — 1 < n. A similar set of 
diagrams can be drawn for the red-shifted components. 
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FIG. 5: The dominant channels for cascading the optical comb 
from a bichromatic pump. Every pair of arrows point from a 
mode in which a pair of photons is annihilated to the modes 
in which the photons are created. 

The coefficients £,n can be determined by calculating 
the powers delivered to a particular mode n via all possi- 



ble channels, and adding these powers up. Notice that by 
adding the optical powers, instead of fields, we disregard 
interference between different channels. This approxima- 
tion may be suitable when the number of such channels is 
large. Requesting that the sum should be equal to P*^"\ 
we arrive at the following recursive relation: 



(n-l)/2 

Sn ^ / ^ 4m4n — 2m — li 
m=0 

which yields a sequence £,n— 1, 1, 2, 3, 5, 8, 16, 27, 
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FIG. 6: A three-FSR-spaced comb generated by two 4.75 mW 
pumps is fit by the envelope functions based on a single- 
cliannel excitation (solid straigfit line) and on multiple ex- 
citations channels (line with dots.) 



Equation (jlOp allows us to fit the comb envelope, see 
Fig. [S] The comb shown in this figure is generated by 
two pumps of 4.75 mW each, set three FSRs apart. Also 



shown are the multi-channel envelope ([T0|) and the expo- 
nentially decaying envelope expected for single-channel 
excitation (when all ^„ = 1). We see that the multi- 
channel model approximates the observed comb much 
better than the single-channel one. 

Let us point out that the processes described by Eq. ([3]) 
will also take place in the single-pump four-wave mixing 
process as soon as the threshold is exceeded and the comb 
arises. Moreover, given that Pq <C Pth, processes ^ will 
dominate the further comb dynamics. In particular, they 
may be expected to build up coherence between the comb 
components, which could not be expected from multiple 
independent processes shown in diagram Fig. [Tf a) for var- 
ious N. A high-purity RF beat note reported in |9l.[lllll3| 
indicates the presence of such a coherence. 

To summarize, we have demonstrated a new and highly 
efficient method of generating optical combs with a rela- 
tively low-Q WGM resonator and only a few milliwatt of 
the external CW pump power. The absence of a pump 
threshold in this process has been confirmed. The theo- 
retical predictions for the shape and pump power depen- 
dence of the bichromatically-pumped comb reasonably 
agree with the experiment. We demonstrated the comb 
spacing varying from one to ten FSRs without any sig- 
nificant change in the comb behavior. This suggests that 
the spacing could be made even larger. Variable spacing 
of the comb may be convenient for the spectroscopy ap- 
plications. We also would like to point out an interesting 
possibility of creating Moire comb pattern, by using three 
or more unequally spaced pump frequencies. This may 
enable an interesting approach to optical synthesizers. 
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